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Abstract. Suspensions of self-propelled particles are studied in the framework of two- 

Y^ [ dimensional (2D) Stokesean hydrodynamics. A formula is obtained for the effective 

• i-H ' viscosity of such suspensions in the limit of small concentrations. This formula includes 

^^. the two terms that are found in the 2D version of Einstein's classical result for passive 

<~i ' suspensions. To this, the main result of the paper is added, an additional term due to 

, MH ' self-propulsion which depends on the physical and geometric properties of the active 

suspension. This term explains the experimental observation of a decrease in effective 
viscosity in active suspensions. 
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O ! 1. Introduction 

00 
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Recently there has been considerable interest in understanding the dynamics of systems 

of active interacting biological agents (active systems, for short), such as flocking birds, 
schooling fish, swarming bacteria, etc [U [21 El IH El [6] . Properties of these strongly self- 
organizing dissipative systems are of fundamental interest to nonequilibrium statistical 
dynamics [3 El El [lOl [IS] and to potential technological applications [13]. 

Common swimming (motile) bacteria, such as Bacillus Subtilis, Escherichia coli, 
and many others, are rod-shape microorganisms (length about 5 fim, diameter of the 
order 1 fim), propelled by the input of mechanical energy at the smallest scales by 
the rotation of helical fiagella attached to the cell wall. Suspensions of swimming 
bacteria (active suspensions) are a convenient representative of self-organizing biological 
systems. At relatively high filling fractions of bacteria, they interact mostly through 
hydrodynamic entrainment induced by their swimming with respect to ambient fluid 
[H [H], which are instrumental in the establishment of long-range order. At the same 
time, they are amenable to accurate experimental studies, control, and manipulation 
|14] . In the dense regime, the dynamics of active suspensions is dominated by the 
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multiple body interactions and long-range self-organized coherent structures, such as 
recurring whorls and jets with the spatial scale exceeding the size of individual bacterium 
by an order of magnitude. This, however, makes the analysis rather challenging. 

In the dilute regime, the hydrodynamic interactions between bacteria (or passive 
particles) are often ignored as the interparticle distance exceeds the range of the flows 
resulting from the particle motion. Therefore, in this case it is possible to isolate the 
effect of the particle-fiuid interactions on the effective properties of active suspensions. 
This constitutes a step towards the full understanding and utilization of the novel 
properties of active suspensions. 

The first step was taken in [15] by Einstein, who derived an explicit formula for 
the effective viscosity of a dilute suspension of passive spheres. This formula shows an 
increase of the viscosity over that of the ambient fluid alone. The correction to the 
effective viscosity is to the first order in the volume fraction of the inclusions. Later, 
Batchelor and Greene obtained the second order asymptotic formula, which takes into 
account pairwise interactions between particles [I6]. Jeffrey calculated the effective 
viscosity of a suspension of ellipsoidal particles to first order in [17] . Unlike the spherical 
case, here the value of the viscosity is affected by the distribution of orientations of the 
inclusions, which is assumed to be uniform. 

In this work, Einstein's classical dilute limit result is extended to the case of 
self-propelled disk inclusions in a 2D Stokesean fluid. The choice of two-dimensional 
hydrodynamics is motivated by the tractable nature of calculations involving Green's 
functions in 2D and also by the quasi-two-dimensional thin film geometry of the 
experiment |14j. In particular, the correction to the effective viscosity is explicitly 
calculated as a function of the orientations of the bacteria's flagella, the intensity of the 
their force of self-propulsion, and the volume fraction of the suspension. In the case of 
an absence of self propulsion and, alternatively, the case of a uniform distribution of 
orientations, the result recovers the 2D version of Einstein's result (see, e.g. [18]). 

The model of the bacterium and the corresponding hydrodynamics are described 
in Section [21 Following Batchelor [19] , the effective viscosity of a dilute suspension is 
defined in Section [21 in terms of a suitable background flow and the disturbance flow 
produced by the inclusion of a bacterium in the background flow. In the case of self- 
propulsion the disturbance flow is due both to the passive response of the background 
flow to the inclusion of a particle as well as the response of the stationary fluid to the 
particle's self-propulsion. Because of the linearity of Stokesean hydrodynamics, the two 
components of the disturbance flow can be computed independently using the Green's 
function for the disk. 

The disturbance flow Up due to the locomotion of a single bacterium is calculated 
in Section [31 This flow translates the bacterium. It is also shown that this disturbance 
flow decays at infinity. This result is emphasized since it is a point of well-known 
deficiency of 2D hydrodynamics. Indeed, while the translation of a ball in 3D produces 
no flow at infinity, in the well-known Stokes paradox, the 3D fiow due to an infinite 
cylinder (described by the corresponding force monopole) translating transversely to its 
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axis generates a non-zero flow at infinity (see, e.g., [21] )• Thus prevents the decouphng 
of rods, which is needed for the dilute limit, where it is assumed that the flow due to 
a suspension can be approximated by the sum of solutions due to a single inclusion. 
However, a self-propelled particle is constrained by the viscous drug force opposing the 
propulsion force, thus producing a force dipole which decays at infinity, unlike the force 
monopole in the case of the moving disk. 

For the other disturbance flow, u[^^, due to the bacterium's passive response to 
a background flow, care must be taken to ensure that no flow is produced at infinity. 
Thus, in Section H] such a flow is selected to ensure the disturbance decays at infinity. 
Once it has been done, the effective viscosity is calculated as a function of the flagella's 
orientations relative to the background flow. The time-dependent nature of the effective 
viscosity, due to active alignment of the inclusions to the flow, is also discussed. 
Conclusions are discussed in Section [5l Finally, the details of calculations are discussed 
in the appendices. 

2. Model and its Homogenization 

A 2D dilute suspension of bacteria is modeled as a collection of discs of radius a, each of 
which has an associated point force, representing the flagellum, placed at a distance Xa 
from the disk, as shown in Figure [H The point force is directed radially outward from 
the center of the bacterium and has some orientation angle a, measured from the x-axis. 
The bacteria are distributed throughout an ambient fluid of viscosity t] which takes up 
the entire domain M^. A single bacterium in an infinite fluid moves with respect to the 
fluid with constant velocity driven by the point force c. 

The suspension is assumed to contain sufficiently many bacteria so as to produce 
appreciable changes in the properties of the equivalent homogenized fluid. At the same 
time, the size of the bacteria and the propulsion force are assumed to be sufficiently 
small, so that the disturbance flow produced by bacteria is negligible at inter-bacterial 
distances and hence can be ignored at the locations of other bacteria. Therefore, in the 
following calculations, it is sufficient to consider a single bacterium with orientation a in 
an unbounded volume of fluid. This is analogous to Einstein's assumptions for passive 
suspensions. 

This unbounded volume is an idealization of the microscopic volume element 
surrounding a single one of the particles located within a macroscopic fluid element. 
Because of the decay assumption, this idealization is permissible for the solution of the 
Stokes' equation within that region. The viscous dissipation in the macroscopic fluid 
element is then the sum of the dissipations in all of the subordinate microscopic volumes, 
containing a collection of particles having some prescribed distribution of orientations. 
Analogous to the method of Batchelor in [19] , the effective viscosity t]* of a suspension 
in an ambient fluid of viscosity 1] is defined as the viscosity of an equivalent fluid with 
no inclusions that produces the same energy dissipation in each macroscopic volume 
element. 
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For a suspension of passive particles, the existence of an equivalent Newtonian 
fluid with the scalar effective viscosity rj* is a classical result of homogenization theory 
(see e.g., [20] and references therein). The analysis for active particles is more subtle 
and has not been carried out in a rigorous mathematical context. The assumption is 
made that the effective viscosity can be defined in the same way for active particles. 
Nevertheless, the technique is presented somewhat cryptically in the literature. To 
clarify the main points, the conceptual side of the calculation is presented here in some 



detail; the technical details can be found in Appendix A (also see [I9] for the derivation 
in the case of passive particles). 



c 



Figure 1. Schematic presentation of a bacterium with orientation a. 

Since the effective viscosity is defined in terms of energy dissipation, a background 
flow (which describes the flow throughout the homogenized fluid and on the boundary 
of the suspension) is chosen so that it has a non-zero rate of energy dissipation (i.e., 
requiring a non-zero strain rate). Choosing M^ to represent the microscopic volume with 
coordinates {xj}, i = 1,2 the velocity components of the background flow are Ui = eijXj, 
where e^j is constant and symmetric strain rate tensor but otherwise left to be specified 
later. Adding a single self-propelled disk inclusion VJ, somewhere in the plane produces 
the disturbance flow u'. Provided that u is selected properly, u' vanishes at infinity, 
decoupling the bacteria, as required of the dilute limit. The disturbance flow can be 
computed explicitly and it is done below. The total flow of the suspension is u = u + u' 
and in principle it should be possible now to compare the dissipation rates of u and u 
by integrating over M? and M? — Vb respectively. 

As shown in [19] , the integration over the unbounded volume M^ leads to technical 
difficulties. This motivates the more efficient method of comparing the dissipation rates 
within a bounded domain V which contains sufficiently many bacteria. Figure [2] shows 
such a domain encompassing two bacteria. Inside V are the regions representing the 
bodies of the bacteria, Vb; and, for each, an e-disk Vs^i around the corresponding point 
force. The disturbance flow in the domain V — J2i ^Bi — Yl ^e,i is required to vanish at 
the outer boundary F = dV. This models a flnite container with the background flow 
applied at the boundary and simplifles the computation of the dissipation rate integrals 
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Figure 2. Bounded domain V for which the modified dissipation rate is computed; 
includes the complement of two a-disks (bacteria) and e-disks (neighborhoods of point 
forces) . 



by replacing a volume integral with an integral over the boundaries of the bacteria. 

Calculating the disturbance flows u', however, involves solving a boundary value 
problem in the two-dimensional domain V — Yli ^Bi — Yl '^e,i- To avoid this, the total 
dissipation rate in this domain can be equivalently calculated as the integral of the 
work of forces at the boundary T . By requiring that u' vanish at F = dV , the flow 
of the homogenized fluid and total flow of the suspension agree on the boundary F, 
but the corresponding values of the stress tensors a and a are not 



Ui 



Ui 






necessarily the same. Now, the density of work at the boundary is riiaijejkXk and 
nidijejkXk for the homogenized and total suspension flows respectively. The effective 
viscosity is defined by integrating these quantities over F and setting them equal: 



CTijUiTijdA = / aijUiUjdA. 
Using the Stokes equation V ■ a' 



(1) 



and repeatedly applying the divergence theorem 
(see appendix Appendix A ), the integrals of these densities over F can be reduced to 
integrals over the boundary of the bacterial domain only: F^ + F^ = dVi, + dV^. These 
integrals are independent of F and upon passing to the limit V — > M? they are expressed 
in terms of the disturbance flow u' vanishing at infinity. 

The effective viscosity is then determined from the expression resulting from 
summing over all inclusions: 

2(r/* - T^)e,,e., (2) 

= ^ X] / {^'ij^kUj - 27]u[nk) dA. 

^ I JdVB^udv.^i 

Applying the dilute assumption, the solution to (lA.Sp is replaced by a superposition of 
solutions for the flow due to a single bacterium, derived in the next section. Thus, for 
aligned bacteria, the above integral can be replaced with one over the boundary of one 
bacterium in M^, 



2(?7* - r])eijeij 



(3) 



Effective Viscosity of Dilute Bacterial Suspensions: A Two- Dimensional Model 6 

eik(f) / {cr'ijXkUj - 2r]u[nk) dA, 

where is the volume fraction occupied by bacteria. 

3. Flow due to a single bacterium 

Let Q be the exterior of the unit disk in M? and 5^^ be the Dirac delta evaluated at 
r — -To, where r*o is the location of the point force. In the low Reynolds number limit, 
the flow due to the bacterium obeys the Stokes equation 

■ijAu'p = Vp'p + c6i^g in fi 

V-u'p =OmQ , (4) 

Up = on r 

where c is the force strength and F = dQ. For simplicity of calculations, set c = cei and 
fo = (1 + A)ei. Note that here, u' is set to on F. This is the reference frame in which 
the bacterium is moving-for the purpose of finding the total flow of the suspension, the 
reference frame in which the fluid is at rest at infinity must be used. 

By writing H as the curl of a stream function $ and taking the divergence of (j4]), 
this reduces to the inhomogeneous biharmonic equation 

$ = |i = on F ■ ^ ^ 

an 

The Green's function G for the domain Q with these boundary conditions is known 
and yields the solution 

c dG 

Hx,y) = — -—{x + zy,l + X). (6) 

Explicitly, in polar coordinates, 

Hr,e)= (7) 



cr cos(6') 
+ log 



((l + A)^-l)(r^-l) 
1 + (1 + A)2r2 - 2(1 + X)r cos(e) 
l + A)2 + r2-2(l + A)rcos(0) 



1 + (1 + A)2r2 - 2(1 + \)rcos{9)J_ ' 

It is easy to see that Up — > Uoo = const, as r -^ cxo. Subtracting this flow (moving 
to the reference frame where the water is at rest at oo) and rescaling to a disk of radius 
a, is shown in figure [3] and is given asymptotically by 

M(r, 6) = (8) 

2gg±^i(cOs(0)+COs(3^))+O(;^) 

Now, by taking the divergence of (j4]), one obtains 
Vp ■ n = r]Au ■ ra on F 
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Figure 3. Velocity field produced by a bacterium in 2D geometry. 



The solution to (jH]) is found using complex function theory (see appendix Appendix D ) 
and is the real part of the complex function 

c a^X\2 + X)^z 



mz) 



27r (1 + A)(a(l + A) - z){{l + X)z - af 



(10) 



4. Calculating the effective viscosity 

In order to calculate the effective viscosity, it remains to perform the integration in ([3]). 
The effective viscosity is independent of the background flow, so one is free to choose 
the flow which simplifies the necessary calculations. In ([3]), u' and p' are the solutions 
to 



r]Au' 
V -u' 



Vp' + C(5r=-o 

= 



subject to the boundary conditions 






^ijXj 



as r — i> 

xeT. 



oo 



:iii 



(12) 



W = const 

Note that, as mentioned in section [31 one is not free to choose the value of the constant 
that H takes on F. 

The formula for the effective viscosity does not depend on the choice of background 
flow (i.e., the choice of the strain rate tensor e^j in (fT2|) ). However, the choice of this 
flow is not completely arbitrary. First, since the effective viscosity is defined in terms 
of energy dissipation, a non-zero rate of strain is required so that there is energy being 
dissipated due to viscosity. Second, the background flow is chosen so that it does not 
rotate bacteria, which allows for the computation of the effective viscosity for a fixed 
(time-independent) distribution of orientations of bacteria (see figure [1]). Third, upon 
placing a disk (i.e., bacterium) in the flow, it is necessary that the disturbance flow 
"^int — > as r ^ cxD. If this were not the case, the dilute assumption that this flow is 
negligible at the locations of other bacteria would be violated. In 2D, this is only possible 
for a background flow that satisfies the conditions for existence of such a decaying flow 
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derived in [22]. The background flow is chosen to be the simplest flow that matches 
these criteria-the sum of two perpendicular shearing flows, 

1 

1 



2en 



(13) 



where eo = const is the amplitude of the strain rate. 

The disturbance flow due to the presence of a bacterium will be 



u 



Up{a) + Ui^t 



with the pressure given by 

p' = p'p{a) + Pint 
where, using notation, u'p{a) := 



u'p{r, 9 



a. 



as a polar vector and p' (a) 



Pr 



(14) 

(15) 
(r, 6'-a), 



rotated to allow for an arbitrary orientation of a bacterium in the background flow. 

For a suspension of bacteria that are not fully aligned, the effective viscosity can be 
written in terms of the function X{a) which gives the probability density of orientations 
of the bacteria a^, which are assumed to be independent identically distributed random 
variables. In this case, ([3]) becomes 

2(?7* - r])eijeij = (16) 

eik(f) / / {a[j{a)xknj 

Jo JdVBUdVe 

—2riu^{a)nk) dAX{a)da, 

where is the volume fraction occupied by bacteria. Performing the integration over 
OVb U dVe, one obtains the following formula for the effective viscosity: 

7]* = r][l + 2(j) (17) 

c 1 + 2A(2 + A) 



+ 



2lT 



sin (2a) X{a)da 



27Tr]eo (1 + A)3 j^ 
where is the two-dimensional volume fraction. 

Note that for isotropically oriented bacteria (X(a) = l/27r = const), the angular 
term goes to zero and the effective viscosity reduces to the result in [18] for the effective 
viscosity of a two-dimensional dilute suspension of disks, 

7]* = r] {1 + 2(f)) . (18) 

Nevertheless, in experiments a reduction in viscosity is anticipated [22]. In reality, 
bacteria are not spherical. An interaction with the background flow thus produces a 
torque on the bacteria which will tend to align them in particular preferred directions 
(see e.g. [21]). For a fully aligned suspension of bacteria, one obtains, from (HT 

c 1 + 2A(2 + A) 



7] = 1] 



+ 



27rr/eo (1 + A)= 



I sin 2a 



(19) 



Stt 



For the background flow flT3|) . the corresponding preferred directions are a = — f , 4 
Thus, after aligning. 



7] = 1] 



c 1 + 2A(2 + A) 

2TT7]eo (1 + A)3 



(20) 
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This formula can be used to explain the experimental data on the reduction of viscosity 
in bacterial suspensions observed in pjj. 

5. Conclusion 

Here a rather puzzling result has been obtained: the effective viscosity of active bacterial 
suspension may decrease with the increase in volume fraction of particles. Moreover, 
the viscosity may even become formally negative (when the flagella add more energy 
than is being dissipated) if the volume fraction and the magnitude of the propulsion 
force c exceed critical values. Manifestations of negative viscosity have been observed 
in experiments [T^ via large-scale instability and the formation of non-decaying whorls 
and jets of collective locomotion (see also simulations [26], [27] for dense suspensions). 
The reduction of effective viscosity can be interpreted as a result of transformation 
by swimming bacteria of chemical energy of the surrounding nutrient medium into 
mechanical energy of fluid motion, and thus replacing energy loss due to viscous 
dissipation. 

It is also notable that the coefficient eg, which characterizes the strength of the 
background flow, shows up in the formula for effective viscosity (unlike in a Newtonian 
fluid). Thus, the additional term due to self-propulsion is inherently non-Newtonian. 
Since eo appears in the denominator, there is a blow up in the effective viscosity as 
eo — *■ 0. Nevertheless, this is to be expected, as the viscosity is, roughly, the ratio of 
stress of the entire flow to rate of strain of the background fiow-in the case of a bacterial 
suspension, as the rate of strain of the background flow goes to zero, the stress does 
not, since the bacteria will continue swimming. 

We also comment that the limit of vanishing background flow eo ^ is not 
trivial. As mentioned above, the reduction of viscosity occurs only when the bacteria 
become aligned in a certain direction determined by the principal axis of the strain- 
rate tensor. However, when eo — > 0, the bacteria no longer have a tendency to 
align (also, a preexisting ahgned state would be destroyed by fluctuations, e.g. due 
to random tumbling of bacteria) and hence the third term in Eq. (ITTll does not 
blow up. Thus, in the limit of eo — > 0, the steady state is in fact the isotropically 
oriented suspensions X{a) = l/27r, which recovers the classical result for disks Eq. 
( TT8|) . Moreover, since the effective viscosity for small strain rates is determined by 
the interplay between the response of the bacterium, through its orientation, to the 
shear strain and random fluctuations, one needs to solve self-consistently the equation 
for the orientation distribution function X{a) in the presence of non-zero strain rate. 
This isotropic orientation can also explain the results of recent simulations [28] in 
which no reduction of the effective viscosity due to swimming bacteria (modeled by 
spherically symmetric particles) was detected: the reduction of viscosity requires particle 
reorientation by shear flow. However, the orientation of spherical particles is not affected 
by shear flow. In contrast, the effect is expected for elongated particles (ellipsoids or 
cylinders with large aspect ratios (of the order of 1:5 for swimming bacteria). 
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Obviously, more dedicated controlled experiments with bacterial suspensions and 
further generalizations of the obtained result to more general flow geometries, such 
as three-dimensional films and slabs, are keenly needed. The puzzling phenomenon 
of reduced viscosity in bacteria-laden fluids may find rather unexpected technological 
applications in bio-medical research and chemical technology, such as microscopic 
bacterial mixers and chemical reactors [T3] . 
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Appendix A. Calculation of Dissipation Rate and Definition of Effective 
Viscosity 

Following [19], the rate of work being done at the boundary T = dQ is given by 

/ aijUiUjdA = Cik / {-PSij + 2r]eij)xknjdA, (A.l) 

where ctjj is the stress tensor, u is the velocity of the fluid, n is the unit outward normal 
for the surface F, P is the pressure of the fluid, and eij is the rate of strain tensor. 
Henceforth, primed quantities will denote the disturbance values due to the presence 
of a bacterial suspension. The new stress tensor can be expressed as dij = aij + a[y 
Additionally, on F, Cij = eij. Thus, the effective viscosity of the suspension i]* is defined 
by setting 

eik / {-PSij + 2r]*eij)xknjdA = (A.2) 

Cifc / {-PSij + 2r]eij + a-j)xknjdA. 

Noting that the terms involving P are identical on both sides and employing the 
divergence theorem yields 

2A{ri* - r])eijeij = tik I a-jXkrijdA, (A. 3) 



r 



where A is the area of the surface F. The right hand side of flA.3P is the additional rate 



of dissipation due to the suspension. Employing the divergence theorem yet again, this 
integral is transformed into an integral over the surfaces of the particles, producing 



eik / cr'-jXkTijdA 



(A.4) 



da',, 

ik 



e^u I ( -^Xk + a[, 1 dV 

+ eifcV' / (T[^XknjdA, 

JdVuyJdVe 
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where Vb is the volume occupied by a single bacterium, V^ is a ball of radius e around 
its corresponding point force, and the summation is taken over all bacteria inside Q. 
Now, the fluid in fl obeys the inhomogeneous Stokes equation 

v-u' =0 mn-j:,VB, ^ 

u' = vi on OVbi 

u' = eijXj on dfl 

where the subscript / has been added to all quantities that can vary among the bacteria. 
In particular, vi is its velocity of the Ith bacterium, r; is the location of its point force, 
and Q indicates the strength and orientation of each point force. Thus, -q^ = in 
^ - E ^B U K. Additionally, 

e^k [ cj[,dV (A.6) 

/■ du'- 



-CifeV" / 2r]u[nkdA, 

JdVnVJdV, 



since the surface integral over F vanishes. Applying these observations to ( ]A.3f) produces 

2{r]* - r])eijeij (A.7) 

= ^ 5Z / i^ij^kUj - 2r]u[nk) dA. 

Appendix B. Velocity 

Recall that the velocity field u solves the Stokes equation, 

riAu = Vp + cSifg in ^2 

V-u =Omn . (B.l) 

u = on r 

For simplicity, it is assumed that c = ex and tq = (A + l)a;. Since V ■ m = 0, the 
velocity can be expressed as the (2D) curl of a scalar stream function $(x, y). This curl, 
which operates on scalar functions and produces a vector function, is defined by 




Curlv^{x,y) 

The scalar curl Curls, which operates on vectors, is defined as 

Curlsu 



dux duy 



dy dx 
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so that CurlsCurlv^{x,y) = A$(x, y). Substituting u = Curly(^{x,y) into (jl]) and 
taking the scalar curl of both sides yields the inhomogeneous biharmonic equation 

A2$(x, y) = Curls-x6r-, = ^%o (B.2) 

Tj oyrj 

with boundary conditions 

(9$ 
$ = ^— = on r. (B.3) 

on 

The Green's function for the domain Q with these boundary conditions is derived 

in the exact same fashion as that for the unit disc, and has a simple form as a function 

of the complex variables z = x + iy and C = ^ + ^7; taken from [25] : 






1-C^ 



(B.4) 



+ li^(N^-i)(lcr-i)- 



This yields the solution formula 

f do 

^ix,y)= / G{x + iy,^ + i-f)- — 5f^^{^,-f)d^d-f 

c f d 

= / -^Gix + iy,^ + i'y)6f^,,d^d'y 

V Jn9-f 

and hence 

c dG 

<!>{x,y) = — -—{x + zy,l + X). (B.5) 

r] 0*7 

Appendix C. Green's Function 

This derivation follows that in [25], with the only difference being that, in this case, the 
domain is the outside of the unit disk D. Nevertheless, the mathematical details are 
identical. The fundamental solution of the biharmonic equation is 

A(^,C) = ^k-CI'iog|;.-CI, (C.l) 

and the general solution of the biharmonic equation is 

T{z) = 2Re{z<!>{z) + '${z)}, (C.2) 

where $ and \E' are arbitrary analytic functions. Thus, the problem is to find a function 
of the form flC.2p that cancels fIC.ip and its normal derivative on dD. To facilitate this, 
flC.2l) can be equivalently written as 

T{z) = 2Re{[z^ -imz) + '^{z)}. (C.3) 

Let T{z, C) = A + T. Then, the condition F = on dD is equivalent to 

Re{v^(^)} = -^(^-C)(^-C)logk-CI (C.4) 

= Re |-i-(z - C) (^ - C) log(l - C^) 
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since z = - there. Since the equation inside the braces is analytic in the unit disk, it 
must be '^(z). Additionally, the condition that F = |^ = on dD implies 1^ = 1^ = 

^ ' -^ ' on ^ ax ay 



there, and so 



IGtt 
Substituting \l/ yields 

Re{$(z)} = 



^ = -zRemz)} + ^-^'{z) + Y^(^ - 
+ -^(^-C)log(^- 0(^-0 = 0. 



(C.5) 



(C.6) 



Re 



t 167r 
and hence, using z = - on dD once more 



^(l--] log(l - C^) 

Stt V z 



$(z) 



ICP 



IGtt 



l(l-f|log(l-C--) 



Thus, 



r(^,C) = — k-crio? 

OTT 



^-c 



c^ 



(C.7) 



(C.8) 



+ 



IGtt 



(kr-i)(icp 



1). 



Appendix D. Pressure 

By taking the divergence of (j4]), one obtains 



Ap 



—V ■ cx6 



r-o 



-c—5^ 



fD.ll 



Vp ■ n = r]Au ■ n on F 

Thus p{x,y) is the real part of some function Il{z), which is holomorphic in ^2^, 
such that n'(z) = T]{Aux — iAuy) on F. In fact, viewing u^ and Uy as functions of x and 

y, 

U{z) = r^j[Au,iX{C),Y{C)) (D.2) 

-tAuy{xh,Y{0)]dC, 

where X{z) = ^ (-2 + ^), ^(-z) = ii (^ ~ 2)' ^^"-^ ^ ^^ ^^^ path from the origin to z. It 
is clear that p = Re (11) satisfies the compatibility condition and is harmonic, except 
where the integrand has singularities, so it remains to check that ARen(z) = —0-^6^^. 
Performing the integration in (1D.2P gives 

-c a^X\2 + X)^z 



^^^^ ' 27r (1 + A)(a(l + A) - z)((l + X)z - a)^ ' 
Taking the real part of flD.3l) and expanding it about x = a{l + A) yields 
, ^, cos(6') ^,_,, 



thus, indeed, Ap{r, 



27rr 



(D.3) 



(D.4) 
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Appendix E. Disturbance flow 

In terms of the stream function of the background flow, '^o{r,9) = r^(cos^6' — sin^6'), 
the stream function for the disturbance flow Ud, from [22] , is given by 



r^ — Ir^a} ^ (a} ^ 

^= 2 ^0 — ,^ 

a* \T 

2 



(E.l) 



a 



+ -Ar^--')l^ol-,e 



,2\2 



4a4 
This yields the velocity 



(r2 - a^ ^2 



r 



r'^o{-,0 



Ud 



r ae ' 
ar 



The corresponding pressure pd is, once more, the real part of 



U{z) = 7^l[AudAX{0,Y{0) 
-tAud,yiXiC),YiO)]dC. 



(E.2) 



(E.3) 



Performing this integration gives 



Pd{x,y) 



IQa^rjxy 



(E.4) 
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